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, ■ Abstract 

^ ' This paper presents a Prolog interface to the MiniSat satisfiability solver. Logic program- 

^ . ming with satisfiability combines the strengths of the two paradigms: logic programming 

' for encoding search problems into satisfiability on the one hand and efficient SAT solving 

I on the other. This synergy between these two exposes a programming paradigm which 

■ we propose here as a logic programming pearl. To illustrate logic programming with SAT 

fSJ ' solving we give an example Prolog program which solves instances of Partial MAXSAT. 

, To appear in Theory and Practice of Logic Programming (TPLP) 

l> ■ 

O ; 

c/5 . 1 Introduction 

O ■ 

^ ' The use of SAT solvers in a wide range of applications is a great success of recent 

• ^ . years and quickly growing more popular. Contributing to this success are two main 

' factors: (a) SAT solvers are continuously becoming more powerful, and (b) encoding 

' . . . 

' techniques are emerging to represent a wide range of search problems as proposi- 

tional formulae such that each satisfying assignment of the encoding represents a 
solution of the problem. 

This paper presents a Prolog interface to the MiniSat SAT solver (Een and 
Sorensson 2004; MiniSAT 2006). MiniSat is a smaU (« 1200 lines of C-code) and 
efficient open-source SAT solver which is designed to enable easy integration with 
other tools and languages. Application of this interface facilitates the best of both 
worlds: Prolog programming for encoding a search problem into a propositional 
formula on the one hand, and application of powerful SAT solving tools on the 
other. 

In recent work (Codish et al. 2006b) we apply SAT encodings to decide LPO 
termination (Dershowitz 1982) of a given term rewrite system t. Here LPO termi- 
nation of r is encoded as a propositional formula ipr and any assignment satisfying 
ip-r indicates a precedence on the symbols in r such that its corresponding lexico- 
graphic path order orients the rules of r. Namely for each rule £ ^ r in t, £ >-ipo r 
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holds. Prolog is well suited for the tasks of parsing the given system t and construct- 
ing the encoding (pr . Searching for an assignment which satisfies (fr is conveniently 
delegated to the SAT solver. 

In other recent works Codish et al. apply SAT encodings to determine termina- 
tion of rewrite systems using dependency pairs (Arts and Giesl 2000) and using 
recursive path orders (RPO) (Dershowitz 1982). A contribution of those works is 
the prepositional encodings for so called argument filterings and usable rules, in 
(Codish et al. 2006), and the encodings for multi-set orders and lexicographic path 
orders modulo permutation of arguments, in (Annov et al. 2006). 

All of the above mentioned results have been developed on top of the described 
Prolog MiniSat interface. Wc have also had positive experience in using the interface 
as an educational tool for teaching logic programming concepts as well as the basics 
of SAT encodings and SAT solving. 



2 Preliminaries 

Most SAT solvers assume as input a prepositional formula in conjunctive normal 
form (CNF). That is a conjunction of disjunctions of literals, or cquivalcntly a 
conjunction of clauses. Each literal is a prepositional variable p or its negation ^p. 
A truth assignment is a mapping from prepositional variables into {0, 1}. 

Syntax: In the logic programming setting we represent: literals as terms of the 
form X or -X where X is a logic variable; clauses as lists of literals; and conjunctions 
of clauses as lists of clauses. For prepositional formulae we use terms involving the 
symbols 0/0, 1/0, -/I, */2, +/2, ==/2, xor/2 for: false, true, negation, conjunction, 
disjimction, bi-implication and xor. The syntax supports also: (X->Y;Z) which is 
equivalent to X*Y+(-X)*Z. In the logic programming setting a truth assignment 
is a substitution of the variables to the constants {0, 1}. We also use a list of 
propositions to represent non-negative integers, in two different ways. The unary 
representation is given by the sum of the bits. For example [0,1,0,1] represents 2. 
The binary representation is as usual, but giving the least significant bit first. For 
example [0,1,0,1] represents x 2° -|- 1 x 2^ -h x 2^ -h 1 x 2^ = 10. 

We now recall the satisfiability and maximum satisfiability problems. Partial 
MAXSAT is a generahzation of SAT and MAXSAT introduced by Miyazaki and 
Iwama (Miyazaki et al. 1996). See also (Cha et al. 1997). 

SAT: Given a prepositional logic formula tp in conjunctive normal form, is there 
an assignment of truth values to the prepositional variables that makes y true. 

MAXSAT: Given a prepositional logic formula ip in conjunctive normal form, 
find an assignment of the logical variables that maximizes the number of clauses in 
if that are true. 
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Partial MAXSAT: Given prepositional logic formulae <^ and ip in conjunctive 
normal form, the problem is to find an assignment that satisfies all clauses in (p and 
maximizes the number of clauses in V that are true. 

3 Interfacing Prolog with MiniSat 

The library implementation is written primarily in SWI-Prolog (Wielemaker 2003; 
SWI-Prolog ) and interfaces the MiniSat solver (Een and Sorensson 2004; Mini- 
SAT 2006) for solving SAT instances. We have integrated MiniSat and SWI-Prolog 
through Ril90 lines of C-code and «140 lines of Prolog code. The C-code handles 
the low-level interface and conversion between the internal data representations of 
SWI-Prolog and of MiniSat. The Prolog code provides a high-level interface for us- 
ing the SAT solver in Prolog applications. The SAT solver is deterministic. It does 
not, and is not intended to, provide alternative satisfying assignments over back- 
tracking. The MiniSat library for Prolog can be downloaded from (Codish et al. 
2006a). It includes three modules for the user: cnf .pi, adder.pl and minisat .pi. 

The module cnf.pl exports the predicate cnf/2. A call to cnf (F, Cnf) trans- 
forms a prepositional formula F to a conjunctive normal form Cnf. The implemen- 
tation of cnf (F.Cnf ) applies a Tseitin transformation (Tseitin 1968) to ensure that 
the size of the conjunctive normal form Cnf is linear in that of the input formula 
F. The basic idea is to associate fresh variables with the sub-formulae of F. So, for 
example, if F is of the form G+H and variables A , B , C are associated respectively 
with sub-formulae F, G, H, then the clauses for the bi-implication A ^ (B V C) 
are introduced to the Cnf together with the clauses for the transformations of B 
^ G and C ^ H. This process establishes the conjunctive normal form of A ^ 
F. Adding the singleton clause [A] then provides a CNF equisatisfiable to F. The 
variables A , B , C are sometimes referred to as "Tseitin variables" . 

A technique proposed by Plaisted and Greenbaum (Plaistcd and Greenbaum 
1986) is applied to reduce the size of the transformation. If it can be determined 
that we are only interested in the truth of a Tseitin variable, we can simplify the 
transformation only including the implication. In the above example, since we are 
only interested in the truth of Tseitin variable A (for F) , the CNF created for A ^ 
(B V C), will be [[-A,B,C]] rather than [[-A,B,C] , [A,-B] , [A,-C]]. Similarly 
when we are only interested in falsity of a Tseitin variable. 

To demonstrate the use of the module consider the following queries where T, Tl 
and T2 are the Tseitin variables. 

?- cnf (X==Y,Cnf) . 

Cnf = [[T], [-X, Y, -T], [X, -Y, -T] ] 
?- cnf ((X*Y)+(-X*Z) ,Cnf) . 
Cnf = [[T], [-T, Tl, T2], [-T2, -X], 
[-T2, Z], [-T1, X], [-T1, Y]] 

Figure 1 illustrates the module adder.pl which includes a textbook (Gormen 
et al. 1990) ripple-carry circuit for binary addition. The module exports the predi- 
cate sum/3. A call to siim(Unary, Binary, Psi) constructs a Boolean circuit Psi. 
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The arguments Unary and Binary are lists of Boolean formulae, the truth values 
of which encode unary and binary numbers respectively. These are the inputs and 
outputs of the circuit which can be seen as computing the binary sum of its unary in- 
puts. The formula, Psi is a propositional statement capturing the correspondence 
between unary and binary representations of non-negative integers. The formula 
Psi is satisfiable exactly when Unary and Binary encode the same number. In the 
recursive call of the code, the formula Psi=Fl*F2*F3 is constructed by splitting the 
inputs into two equal halves (padding with a zero if necessary). From the recursive 
calls, the subformulae Fl and F2 relate the two halves of the inputs (unary) to 
the (binary) numbers SI and S2. The call to add(Sl ,S2,Sum,F3) constructs the 
formula F3 for the binary Sum of SI and S2. Observe that the length of Binary is 
[log2 I Unary |] . 

:- module (adder, [siun/3] ) . 

sum(+,-,-). 
sum([B] . [S] .(S==B)) . 
sum( [Bl,B2|Bs] ,Sum,Fl*F2*F3) :- 

split ( [Bl , B2 I Bs] , Xs , Ys) , 

suin(Xs,Sl,Fl) , sum(Ys,S2,F2) , add(Sl,S2,Sum,F3) . 

split (+,-,-) . 
splitC [],[],[]). 

split([X],[X],[0]). 

split([X,Y|XYs] , [XlXs] , [YlYs]) :- split (XYs, Xs, Ys ) . 
add(+, + ,-,-) . 

add( [X I Xs] , [Y I Ys] , [Z I Zs] , (CXY==CarryXY) * (Z==SiamXY) *Sum) : - 
half adder (X,Y,SumXY,CarryXY) , 
adder (Xs , Ys , CXY , Zs , Sum) . 

adder(+, + ,-,-) . 
adder ([] , [] , Carry, [Z] ,Z==Carry) . 

adder([X|Xs] , [Y|Ys] , Carry, [Z|Zs] , (CXY==CarryXY)*(Z==SumXY)*Rest) :- 
f ulladder (X , Y , Carry , SumXY , Carry XY) , 
adder (Xs,Ys, CXY, Zs, Rest) . 

y. f ulladder (+, + , + ,-,-) . 

fulladder(X, Y, C, (X xor Y xor C), (C -> X+Y ; X*Y) ). 

half adder (+, + ,-,-) . 
halfadder(X, Y, (X xor Y) , X*Y ). 



Fig. 1. A ripple-carry adder circuit 

To demonstrate the use of the module consider the following query construct- 
ing the circuit Psi for the inputs Unary = [X+Y , X*Y , X==Y , X xor Y] and outputs 
Binary. 

?- suin([X+Y,X*Y,X==Y,X xor Y] .Binary, Psi) . 
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Binary = [SI, S2, S3] 

Psi = (T1==X+Y)*(T2==(X==Y))*(T3==T1 xor T2)*(T4==T1*T2)*(T5==X*Y)* 
(T6==X xor Y)*(T7==T5 xor T6)* (T8==T5*T6)*(S1==T3 xor T7)* 
(S2==T4 xor T8 xor (T3*T7))* (S3==(T3*T7->T4+T8;T4*T8) ) 

Table 1 illustrates the declarative meaning of the predicate depicting the truth 
values for Unary and Binary determined by Psi for the 4 truth assignments of X 
and Y. 



X 


Y 


Unary 


Binary 








[0, 0, 1, OJ 


[1, 0, OJ 





1 


[1, 0, 0, 1] 


[0, 1, 0] 


1 





[1, 0, 0, 1] 


[0, 1, 0] 


1 


1 


[1, 1, 1, 0] 


[1, 1, 0] 



Table 1. Truth values for the circuit summing Unary =[X+Y,X*Y,X==Y,X xor Y] 

The module minisat.pl exports four predicates: 

• solve (Cnf) succeeds if and only if the formula Cnf in conjunctive normal 
form is satisfiable, binding its variables to truth values (false) and 1 (true) 
that satisfy Cnf. 

• sat (Cnf ) succeeds if and only if the formula Cnf in conjunctive normal form 
is satisfiable. It is similar to solve (Cnf ) but does not bind any variables. 

• minimize (Vec, Cnf ) is similar to sat (Cnf). The additional argument Vec is 
a list of variables (e.g., occurring in Cnf). The variables of Vec are assigned 
the truth value that minimizes the binary number represented by Vec for all 
solutions of Cnf. If Cnf has no solutions (i.e. is unsatisfiable) it fails. 

• maximize (Vec, Cnf ) is similar to minimize (Vec, Cnf ) but the assignment 
returned maximizes the value of the non- negative integer represented by Vec. 

The predicates minimize/2 and maximize/2 are illustrated in Figure 2. Consider 
the query ?- maximize (Xs , Cnf ) where Xs is a list of k variables and Cnf a formula 
in conjunctive normal form. To solve the query, the basic idea is to pose k questions, 
one for each variable of Xs, to the SAT solver to determine the maximum value 
of Xs. Each question determines the satisfiability of Cnf when setting the next bit 
in Xs to its maximal value. Observe, in the code, that these questions are posed 
using sat so as not to bind variables in the formula. Note also that each of these 
questions is of size 0(|Cnf |). 

To demonstrate the use of the module consider the following queries: 

• solve/ 1. The call succeeds binding X and Y to truth values. 

?- cnf (X==Y,Cnf) , solve(Cnf). 
X=0, Y=0 

• sat/1. The call succeeds without binding X and Y. 

?- cnf (X==Y,Cnf) , sat(Cnf). 
Yes 
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% minimize (+,+). 

minimize ( [] , CNF) : - sat (CNF) . 

minimize ( [B IBs] , CNF) :- minimize (Bs, CNF) , ( (B=0, sat (CNF)) -> true ; B=l ). 

% maximize (+,+) . 

maximize ( [] , CNF) : - sat (CNF) . 

maximize ([B I Bs] , CNF) :- maximize (Bs , CNF) , ( (B=l, sat (CNF)) -> true ; B=0 ). 



Fig. 2. Minimization and maximization in minisat.pl 

• suin/3 with solve/ 1. The first call succeeds binding X=0 and Y=0. In this case 
the circuit output is 1 (only one of the inputs is true under this assignment). 
The second call indicates that it is possible to satisfy 2 of the inputs under 
the assignment X=0 and Y=l. 

?- sum([X+Y,X*Y,X==Y,X xor Y],Suin,F), cnf(F,Cnf), solve(Cnf). 
X = 0, Y = 
Sum = [1, 0, 0] 

?- sum([X+Y,X*Y,X==Y,X xor Y] , [0, 1 , 0] ,F) , cnf(F,Cnf), solve(Cnf). 
X = 0, Y = 1 

• sum/3 with maximize/2. The answer Sum= [1,1,0] indicates that it is possible 
to satisfy at most throe of the four formulae. The call maximize (Sum, Cnf) 
binds only the elements of Sum. To obtain the maximizing assignment, X=l, 
Y=l, the call to maximize must be followed by a call to solve(Cnf ). 

?- sum([X+Y,X*Y,X==Y,X xor Y],Sum,F), cnf(F,Cnf), 
maximize (Sum, Cnf ) , solve(Cnf). 
Sum= [1,1,0] 
X=l, Y=l 



4 Solving Partial MAX-SAT 

To solve a Partial MAXSAT problem for CNF formula ip and we seek an assign- 
ment that satisfies and maximizes the number of clauses of ijj which are satisfied. 
We solve the more general problem with (p an arbitrary prepositional formula and 
ip a list of n propositional formulae. 

The solution is illustrated in Figure 3. The arguments Phi and Psi correspond to 
(p and i/j respectively. The key idea is to construct the formula SumPsi representing 
the sum of the n formulae in Psi. This results in a vector Max with [log2n] bits. 
We then aim to satisfy the conjunction Phi*SumPsi while maximizing the number 
represented by Max. Solving the query maximize (Max, Cnf ) involves O(logn) calls 
to the SAT solver. 

For example the following query provides an assignment which satisfies ip = X+Y 
and 4 of the 7 formula in the second argument. 
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7. partialMaxSat(+,+) . 
partialMaxSat(Phi,Psi) :- 

sum(Psi,Max,S\imPsi) , cnf (Phi*SumPsi ,Cnf ) , 

maximize (Mcix, Cnf ) , solveCCnf ) . 

Fig. 3. Solving Partial MAXSAT 

?- partialMaxSat(X+Y, [X*Y,X==Y,X xor Y,-X+Y, -X, -Y, X]). 
X = 1, Y = 1 

5 Conclusion 

We define a Prolog library for solving SAT instances through an interface to the 
MiniSat solver. The combination of Prolog for manipulating Boolean formulae, and 
powerful SAT tools for solving them is compelling. Wc can straightforwardly build 
solutions to diSicult problems, with small and clean, pearlish code. We illustrate 
this in the paper by encoding a reduction from Partial MAXSAT to SAT. We have 
successfully used the library to solve many such problems (Codish et al. 2006b; 
Codish et al. 2006; Annov et al. 2006). 
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